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Chemical abundance patterns - fingerprints of 
nucleosynthesis in the first stars 
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Abstract. The interstellar medium of low-metallicity systems undergoing star formation will show chemical abundance in- 
homogeneities due to supernova events enriching the medium on a local scale. If the star formation time-scale is shorter than 
the time-scale of mixing of the interstellar matter, the inhomogeneities are reflected in the surface abundances of low-mass 
stars and thereby detailed information on the nucleosynthesis in the first generations of supemovae is preserved. Characteristic 
patterns and substructures are therefore expected to be found, apart from the large scatter behaviour, in the distributions of stars 
when displayed in diagrams relating different element abundance ratios. These patterns emerge from specific variations with 
progenitor stellar mass of the supernova yields and it is demonstrated that the patterns are insensitive to the initial mass function 
(IMF) even though the relative density of stars within the patterns may vary. An analytical theory of the formation of patterns 
is presented and it is shown that from a statistical point of view the abundance ratios can trace the different nucleosynthesis 
sites even when mixing of the interstellar medium occurs. Using these results, it should be possible to empirically determine 
supernova yields from the information on relative abundance ratios of a large, homogeneous sample of extremely metal-poor 
Galactic halo stars. 
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1. Introduction 

During the last decade, the methods for abundance determina- 
tions for faint Pop. II stars have reached a stage when accurate 
abundance ratios can be obtained for great samples of stars with 
overall metallicities below 1/100, and even 1/1000, of the so- 
lar. This has made it possible to explore not only abundance 
trends, e.g. the variation of oxygen or magnesium abundances 
with decreasing iron abundance, but also the intrinsic scatter in 
abundance ratios such as Eu/Fe, at a given Fe/H as a function 
ofFe/H. 

An early example of a detailed discussion of the disper- 
sion in abundance ratios was that of Edvardsson et al. (1993) 
in their study of the chemical evolution of the Galactic disk. 
These authors found a significant scatter in Fe/H for Disk stars 
of a given age and a given galactocentric mean distance; how- 
ever, they did not find any tendency for e.g. Mg/Fe or other 
a-element/iron abundance ratios to scatter at a given overall 
metallicity, except for the most metal-poor stars. A tendency 
for these latter stars could be interpreted as the result of a 
greater star-formation rate (SFR) in the inner Galaxy. In their 
study of Pop. II dwarfs Nissen at al. (1994) found that the scat- 
ter in abundance ratios Mg/Fe, Ca/Fe or Ti/Fe for their sample 
of Galactic halo stars with —3.2 <[Fe/H] < —1.8 was less 
than 0.06 dex. They also found an upper limit in the scatter 
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in O/Fe of 0.15 dex. Since the yields of these different ele- 
ments are different from supemovae (SNe) of different masses, 
the authors could conclude from the small abundance scatter 
that the elements observed in the stars must be the results of at 
least about 20 SN explosions, otherwise statistical fluctuations 
in abundances should have been present. 

The fact that a scatter in abundance ratios for the most 
metal-poor stars should result from a small number of SNe, 
with different masses and therefore different yields of heavy 
elements, was also pointed out by Audouze & Silk (1995). 
McWilliam et al. (1995) explained the large scatter in s-process 
element abundances like Ba/Fe and Sr/Fe for [Fe/H]< —2 
in similar terms (see also McWilliam et al. 1996; McWilliam 
1998; Ryan et al. 1996). 

The abundance scatter for Halo stars has also been mod- 
elled. In a stochastic Halo formation model Argast et al. (2000) 
studied the scatter in relative abundances as resulting from the 
small number of SNe contributing to the abundances for the 
most metal-poor stars; these authors found that a great scatter 
should be expected in several ratios for stars with [Fe/H]< —3, 
representing early evolutionary phases when the interstellar 
medium (ISM) of the Halo was unmixed and dominated by lo- 
cal inhomogeneities. In the range where [Fe/H] increases from 
—3 to —2 there is a gradually diminishing scatter due to the 
contribution from an increasing number of SNe and more mix- 
ing in the ISM, while for still higher metallicities the Halo 
ISM is well mixed, and homogeneous abundance ratios result. 
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Recent studies, empirical as well as theoretical, of chemical in- 
homogeneities in the Galactic halo have also included work on 
r-process elements (notably Eu, as well as Ba, see Ishimaru & 
Wanajo 1999; Raiteri et al. 1999; Travaglio et al. 2001 and ref- 
erences given there) and Be and B (e.g. Suzuki et al. 1999). In 
particular, Tsujimoto et al. (2000 and references therein) have 
developed a stochastic chemical evolution model to study the 
effects of inhomogeneous r- and s-process element enrichment 
in the early Galaxy. Some of their relative abundance diagrams 
show features resembling those to be discussed here. 

There are today three identified metal-poor Halo stars with 
dramatic r-process signatures where CS 22892-052 is the most 
famous one (Sneden et al. 1996). The identification of such ex- 
treme outliers gives important clues to the understanding of the 
nucleosynthesis in these early SN explosions. However, it is 
not possible to determine, from the outliers alone, which type 
of SN (i.e. what progenitor mass) is able to produce such an 
abundance signature, nor to estimate the relative significance 
of different SNe with other progenitor masses. In this paper, 
we shall present an alternative approach to unveil the statisti- 
cal properties of SNe, taking into account a whole population 
of extremely metal-poor stars. The effects from a small num- 
ber of SNe affecting the chemical composition of these stars 
are explored. We shall demonstrate the probable presence of 
fine structure patterns in the diagrams where abundance ratios 
are plotted relative to each other, reflecting the number of SNe 
contributing and the sometimes strong mass dependence of the 
yields. In fact, we shall argue that these patterns, if observed, 
could be used for empirically exploring such properties fur- 
ther. Some preliminary results of our work were published in 
Karlsson & Gustafsson (2000). 

After some general comments and definitions in Sect. ^ 
we shall present our simulations in Sect. ||. The origin of the 
patterns is further analysed in Sect. ^ wherein an analytical 
theory of the distributions of stars in different abundance planes 
is developed. Observational implications are discussed in Sect. 
H and the conclusions are presented in Sect. ^. In the Appendix, 
we derive some general expressions describing the statistics of 
random variables. 

2. Chemical enrichment in metal-poor systems 

Local chemical inhomogeneities in the interstellar medium, 
caused by the first generations of SNe in the Galaxy, may or 
may not be preserved in subsequent generations of stars, de- 
pending on how efficient the mixing of the ISM is relative 
to the rate of star formation. The global mixing can be de- 
fined in terms of a mixing efficiency parameter, emix, which 
describes how fast the mixing occurs as measured in, say, solar 
masses per million years. Hence, Tmix oc Msystcm/emix, where 
Mnystem IS the total mass of the system. The mixing efficiency 
parameter can also be used to define a characteristic mixing if A/mix ^ M. 
mass, such that A/i„ix = emixTsF, where tsf oc SFR^^ is the 
star formation time-scale. 

Thus, if the mixing is efficient enough or/and the star- 
formation rate is low, the time-scale of global mixing is much 
shorter than the star formation time-scale, or equivalently. 
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Fig. 1. a) A sequential build-up of the C/Mg ratio. The ratio 
changes when the contribution from another SN is added. The 
number of SNe is denoted by n. b) The chemical track in the 
[Mg/H]-[C/Mg] plane. The amount of magnesium is mixed 
with 5 X 10^ A4q of hydrogen in calculating the Mg/H ra- 
tio, c) The chemical track in the [Mg/Fe]-[C/Mg] plane, d) 
The variation of the stellar yields with SN mass taken from 
the zero-metallicity models of Woosley & Weaver (1995). The 
pluses (+) denote the carbon yield, the open circles (o) de- 
note the magnesium yield and the squares (□) denote the iron 
yield. Note the downturn in the carbon yield at ~ 25 A^©. At 
the same mass, the Mg-yield increases. This is reflected by the 
jumps in the chemical tracks 



geneities are wiped out before subsequent generations of stars 
are formed and the system will be considered well mixed at 
any instant of time. The stars will reflect the cumulative build- 
up of the elements and the whole stellar population has, in 
this case, a common chemical history meaning that a time-axis 
can be defined. This, in turn, implies the existence of an age- 
metallicity relation. On the other hand, if Tmix ^ tsf, i.e.. 



system 



the stars would trace the local fluctua- 



is comparable to M. 



system- 



In this case the inhomo- 



tions caused by individual enrichment events and the chemical 
inhomogeneous state of the system would be resolved and pre- 
served by low-mass (yet unevolved) stars. Such a system can 
not be considered evolving with time in a unique way since the 
stars in the system describe different chemical histories. Hence, 
a single age-metallicity relation does not exist. 
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Let us discuss in more detail how the chemical enrichment 
of the ISM can be reflected in the stellar population. Suppose 
that we have a system which initially consists of primordial 
gas (i.e. Zinit = 0). Assume that the instantaneous recycling 
approximation is valid and that the whole system is well mixed 
at all times, i.e. Tmix <C tsf- Any massive star that explodes 
immediately pollutes the whole system and the chemical com- 
position of the ISM is changed accordingly. Let us call a chain 
of such enrichment events a chemical series. Note that the pro- 
genitor masses of the SNe do not follow a decreasing sequence, 
they are randomly distributed. A finite part of a series, say, up 
to i SNe, shall be referred to as a chemical sequence which de- 
scribes the chemical state of the system at the time when i SNe 
have exploded. A star formed out of this gas is a realization of 
the chemical sequence. Now, if low-mass stars are continuously 
formed over time the complete enrichment history of the ISM 
is mapped. The chemical evolution of the system is then fol- 
lowed by displaying the stars in different abundance diagrams. 
Such time-lines shall be called chemical tracks and are realiza- 
tions of the chemical series up to a certain number, n, of SNe. 
Since the ISM is homogeneous there exists only one chemical 
series in the system. This is illustrated in Fig. |l]by the sequen- 
tial build-up of the C/Mg ratio. Two chemical tracks related to 
the C/Mg ratio are also shown (Fig. [l|b, [l|c). The distinct jumps 
in Figs. |l]a and |l]b reflect the enrichment of the ISM by a mas- 
sive SN. The yields as functions of progenitor mass are shown 
in Fig. p. 

Systems like the one discussed above show unique chem- 
ical tracks for each pair of elements (or element ratios). This 
is clear since the chemical state of such a system at any time 
is described by a single chemical series. Suppose now that 
-^^inix -A/systcm, 1-6. our systcm consists of a number of 
separate subsystems, or star-forming regions. As above, the in- 
stantaneous recycling approximation is valid for each region 
but the different subsystems evolve differently, and their chem- 
ical states are described by different chemical series. A pop- 
ulation of stars, randomly picked from different star-forming 
regions would, then, describe the chemical evolution of sev- 
eral systems, leading to a large star-to-star scatter in the ob- 
served abundances. However, due to the many different enrich- 
ment histories of these stars, substantially more information on 
the elemental production sites (described as stellar yields) is 
preserved and we shall see that the distribution of stars in the 
abundance diagrams show structures and patterns created by 
specific variations in the SN yields with progenitor mass which 
can provide detailed clues to the production of elements in the 
early Galaxy. 

3. Simulations of the cliemical patterns 

The existence of a large scatter in abundance ratios for the 
Galactic halo stars with [Fe/H]< —2.5 suggests the second sce- 
nario (i.e. Tinix ^ Tsf) discussed in Sect. |2| Below, we shall 
simulate the chemical enrichment of such a metal-poor system. 
In all our simulations particular data on SN yields will be ap- 
plied. These data are rather uncertain, not the least as regards 
the variations of the yields with stellar mass, which also makes 
our actual predicted abundance patterns uncertain. Our point 



here is, however, not to make exact predictions but to illustrate 
the general effects of mass-dependent yields. 

3.1. General assumptions 

The gaseous medium of the system (i.e. the Galactic halo) is 
assumed to be initially primordial, i.e., Zinit = 0. The char- 
acteristic mixing mass is assumed to be much smaller than the 
total mass of the system, such that Mmix Afsystcm- This is 
realized by assuming M,„ix to be of the same order of magni- 
tude as the localized star-forming regions of, say, ^ 10^ Mq 
each. From the number of Halo horizontal branch stars, Binney 
& Merrifield (1998) estimate the present day mass of the stellar 
component of the Halo to ~ 10^ -Mq, which implies that the 
total mass of the Halo was, at least, on the order of Afsystom = 
10^ A4q in the early epochs. Thus, Mmix/A/systcm ~ 10^'^. 

We shall distinguish between abundances measured rela- 
tive to hydrogen, and abundance ratios between different, more 
heavy elements. To trace individual SNe or groups of SNe us- 
ing a ratio of some element A relative to hydrogen (an A/H ra- 
tio), requires knowledge about the mixing of the star-forming 
region after the SNe ejection. Ratios between heavier elements, 
both produced in a SN, are less sensitive to the mixing. This 
is due to the fact that an A/H ratio, in the first-order ex- 
pectation, essentially depends on the mean distance from the 
contributing SNe while ratios between different, heavier ele- 
ments do not. In our models we first schematically assume that 
the newly synthesized elements mix with a constant hydrogen 
mass, i.e. all star-forming regions are supposed to have equal 
mass. Furthermore, the mixing with the cloud material is as- 
sumed to be instantaneous and complete. Subsequently, low- 
mass star formation occurs within each cloud. This implies 
that the time-scale of mixing within the region is shorter than 
the formation time-scale for individual stars which in turn is 
shorter than the epoch of the star formation activity in the re- 
gion and the life-time of the region itself. 

Recurrent SN explosions may disrupt the regions before lo- 
cal mixing and subsequent star-formation occurs. If this hap- 
pens inter-cloud mixing has to be taken into account. This will 
naturally affect abundances relative to hydrogen in the ISM. 
However, abundance ratios between more heavy elements are, 
at least statistically, almost unaffected by such a mixing. 

Subsequently, we shall discuss two types of abundance- 
ratio diagrams (in the general case denoted "A/ diagrams" be- 
low). Some diagrams will relate the abundances of two heavier 
elements to hydrogen, such as [C/Mg] vs. [Mg/H]; these are 
denoted "A/H diagrams". These diagrams will be affected by 
the not very well-known degree of mixing of SN remnants with 
other, less processed, material in the different star-forming re- 
gions. Therefore, we shall primarily discuss diagrams display- 
ing the relative abundances between three heavier elements, 
e.g. [C/Mg] vs. [Mg/Fe], diagrams that are much less sensi- 
tive to the mixing uncertainties. This type of diagrams will be 
denoted "A/ A diagrams" below. 

The ejected material from each individual SN is assumed 
to be well mixed prior to subsequent star formation, such that 
any stars formed out of this material will all have the same 
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abundances. Both observations (see Travaglio et al. 1999) and 
hydrodynamical simulations (e.g. Kifonidis et al. 2000) of 
core collapse supernovae show evidence for substantial mixing 
of elements shortly after the core bounce (however, see also 
Hughes et al. 2000; Douvion et al. 1999 where they discuss 
heterogeneous mixing in the Cas A supernova remnant). The 
simulations also indicate that the interaction between the pro- 
cessed material and the outer layer of hydrogen (i.e in SNe type 
II) causes a complete homogenization of the ejected material. 
This process does not work for stars without a massive hydro- 
gen envelope (i.e. SNe Ib/Ic) for which the mixing might be 
less pronounced (Kifonidis, private communication). 

In our investigation we assume that the stellar yields are 
one-dimensional functions of progenitor mass of the SN. For 
convenience, we also assume that stars less massive than 
10 A^o do not produce heavy elements within the time-scales 
considered here and that the prompt enrichment of the primor- 
dial ISM by a population of very massive stars (> 100 A4q) 
suggested by, e.g., Wasserburg & Qian (2000) did not oc- 
cur. Recent yield calculations by Umeda et al. (2000) show 
a moderate metallicity dependence in the lowest metallicity 
regime for the secondary elements, notably ^^N, ^"^Na and 
^^Al. However, primary elements such as ^^C, ^^O and ^^Mg 
are almost independent on metallicity. A high dependence en- 
ters first after the amount of metals initially present is sufficient 
to cause extensive mass-loss through radiation-driven winds 
(Maeder 1992; Portinari et al. 1998). The chemical yields may 
also be altered by stellar rotation (Heger & Langer 2000). It is 
straight forward to incorporate rotational dependent yields in 
the simulations, as well as in the analytical theory presented 
below (Sect. 4.2). The situation is different for metallicity de- 
pendent yields since the metal content in every star is coupled 
to the history of chemical enrichment. We shall neither con- 
sider metallicity- nor rotational dependent yields in the present 
study. 

For the general discussion we shall assume that the amount 
of any heavy element ejected by a SN of a given progenitor 
mass is constant from one stellar generation to the next. The 
ejecta of an element A, eA{m, Z), can be written as 



still other and different regions. This assumption is further dis- 



cussed in Sect. 3.4.4 and there found to be reasonable. 



eA{m, Z)^{m- nir) x Z™'(< - t„i) + PAim), 



(1) 



where m is the progenitor mass of the star, rrir is the remnant 
mass, (t — Tm) is the initial mass fraction of the element at 
the time of formation of the star, t„i is the life-time of stars with 
mass m and pA {m) is the stellar yield, the mass of element A 
produced or destroyed in the star Hence, the total ejected mass 
is a sum of the initially present amount of A and the newly 
synthesized amount. For low metallicities the first term on the 
RHS in Eq. ^) becomes negligible and, thus, the ejected matter 
is completely dominated by the stellar yield. 

We assume the sampled Halo stars to be statistically in- 
dependent. Thus, we assume that all chemical sequences (i.e. 
stars) are picked randomly from an infinite number of chemi- 
cal series. This assumption means physically that the stars that 
originate from star formation regions where, say, three SNe 
have exploded sample different such regions, and that stars 
coming from regions where four SNe have exploded sample 



3.2. Model I 

To simulate the distribution of abundances in Halo stars we 
adopt the simple picture that star formation and heavy ele- 
ment enrichment of the ISM are confined within star-forming 
regions of about a Jeans mass each, i.e. a hydrogen mass of 
AIyi = 5 X 10^ Mq- This particular value is chosen ad hoc. 
Another hydrogen mass would simply introduce a shift in the 
abundances relative to hydrogen. We place a number of high- 
mass stars with randomly distributed masses (according to the 
Salpeter IMF if nothing else is stated) between 10 and 100 Mq 
in the hydrogen clouds and let them explode as core collapse 
SNe. The total number of SNe exploding in each cloud vary in 
a range from one to a large number (i.e. up to ^ 20). The SNe 
produce heavy elements according to theoretically calculated 
stellar yields. We shall use the yields by Woosley & Weaver 
(1995) calculated from their metal-free models (Z models) and 
by Nomoto et al. (1997). Hereafter, these works will be referred 
to as WW95 and Netal97, respectively. All yields are extrapo- 
lated by a constant to the ends of the mass interval. The yields 
by WW95 are modified according to the possible decay chan- 
nels of the unstable isotopes. The high-mass stars explode di- 
rectly and their ejecta are instantaneously mixed with the cloud 
material. Hence, these regions are chemically different accord- 
ing to the number and masses of SNe that have formed and 
enriched the regions. Subsequently, low-mass stars are formed 
out of the enriched gas. Effectively, this means that we sum 
up the heavy element contribution from every SN within each 
cloud. The abundance of any element A relative to hydrogen 
(by number) is then calculated and taken as the surface abun- 
dance of a low-mass star formed in the cloud. 

We assume the total number of low-mass stars formed in 
each cloud to be constant. This implies that there is an equal 
probability to pick a star tracing one chemical series as it is to 
pick a star tracing another series, and since the number of star- 
forming regions is assumed to be large (i.e. the stars are sta- 
tistically independent) we are allowed to select the stars from 
different regions. The total number of Halo stars (N) in a sim- 
ulated sample is then governed by 



N 



(2) 



1=1 



where Wi is the number of clouds in which i SNe explode and 
n is the maximum number of SNe that is allowed to explode 
and enrich a single cloud. In a typical simulation we assume 
Wi to be constant, e.g. Wi — W — 25 for all i, thus the total 
number of stars is = n x W. By displaying various abun- 
dance ratios for low-mass stars in different A/ diagrams we are 
then able to follow the early chemical enrichment phase of an 
initially metal-free system. 
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Fig. 2. Different mass distribution functions as described by 
Eq. (0). The full line denotes the relative number density of ex- 
ploded stars at < = 28.7 Myrs (corresponding to the life-time of 
a 10 AIq star) for a constant SFR. The dashed line denotes the 
relative number density for an exponential SFR with tsF = 10 
Myrs. The dotted line is the Salpeter IMF 

3.3. Modem 

We have alternatively modified our simple chemical enrich- 
ment model by introducing continuous star formation in the 
regions. This mainly affects two parameters, the mass distribu- 
tion function (which is, in Model I, equal to the IMF) and the 
lower mass limit of this distribution. 

In Model I, the high-mass stars are formed in an initial 
burst. By allowing the stars (both high- and low-mass stars) to 
form continuously over a certain period of time we can account 
for a mass distribution function that changes with time as stars 
with different masses have different life-times. We shall adopt 
an exponential star formation rate (tp) such that 



(3) 



where ^sf is a characteristic time of the star formation. This 
time is set to the same value in all clouds. Furthermore, sup- 
pose that the IMF is a power law of the form (j){m) — 0oto^" 
and introduce r™, the life-time of a star with mass m. The ex- 
pression for the time-dependent mass distribution function (see 
Eq. (7.4) in Pagel 1997) describes the distribution of stars over 
m at a certain t. However, we are more interested in the stars 
that have enriched the cloud at time t. The corresponding mass 
distribution of exploded stars, (see Fig. 0) is governed by 



f{m,t) 



0, 

/om""(l - 



(t-T„)/ts 



> t 



(4) 



where /o — (pQipotsp- The stellar life-times are adopted from 
the lowest metallicity models of Portinari et al. (1998). Note 
that in Model I, the function /(m,t) was equal to the IMF 
since all the high-mass stars were allowed to explode before 
the formation of the low-mass stars. 

The star formation period is assumed to be ~ 30 Myrs in 
each cloud. This is only slightly longer than the estimate by 
Shull & Saken (1995), for OB associations. A lower life-time 
of the clouds would hinder the formation of stars that could 



be enriched in elements produced by the least massive SNe, as 
our models do not take into account global mixing and a second 
generation of star-forming regions. 

Except for the continuous star formation rate Model II is 
based on the same assumptions as Model I. We adopted a con- 
stant star-formation rate in each cloud, which corresponds to 
a long characteristic time, tgp (see Fig. ||), leading to a mass 
distribution of SN progenitors that is significantly different 
from the IMF. The slope of the IMF was set to a = 2.35. 
Furthermore, a read-off time (<rcad-off ), distributed according 
to the SFR, was generated in each cloud which determined the 
actual number of polluting SNe, i, via an integration of the SFR 
(up to troad-off ) normalized to the total number of high-mass 
stars in the cloud. The total number of high-mass stars formed 
in each cloud was randomly generated according to a Gaussian 
distribution centred at n = 10 and with a — 20 (for n < the 
probability is zero). Thus, the number of stars that have been 
enriched by i SNe (Wi) is not a priori known but is determined 
after the simulation stops. As mentioned above, this read-off 
time also sets the lower cut-off of the distribution function of 
the exploded high-mass stars as no star with a longer life-time 
than <rcad-off has been able to enrich the cloud. 

3.4. Results 

We should emphasize that a realistic modelling of the early 
chemical enrichment requires a more physical treatment of the 
mixing than adopted here. The abundances relative to hydro- 
gen are particularly sensitive. The situation is quite different 
for abundance ratios, however They are independent of varia- 
tions in the mixing mass. They are also insensitive to global, 
inter-cloud mixing, infall and subsequent generations of star- 
forming regions (see the discussion in Sect. 43). Therefore, 
we shall, based on this distinction, separately discuss the two 
corresponding types of A/ diagrams. 

3.4.1. A/H diagrams 

Let us first discuss some properties of the A/H diagrams be- 
fore we turn to the A/ A diagrams. We see from Fig. ^ that the 
appearance of the stars in the diagrams depends naturally on 
the produced amount of carbon and magnesium in the massive 
stars. However, the shape of the stellar yield functions are re- 
sponsible for possible trends and/or the groupings of stars into 
different substructures and patterns. These structures appear as 
a result of the various enrichment histories of the low-mass 
stars (cf. Fig. |l| and the discussion in Sect. |^. Even though we 
generally do not expect patterns to appear in A/H diagrams dis- 
playing real observations the overall distribution of stars, such 
as large-scale trends, can still hold important information. If 
patterns would really be observed in A/H diagrams this would 
put strong constraints on the star formation and mixing pro- 
cesses in the early Galaxy. There are three effects that are di- 
rectly observed in the diagrams in Fig. ^ and Fig. Q 

Firstly, the star-to-star scatter seems to decrease with in- 
creasing metallicity (represented by [Mg/H]). This is best seen 
in Fig. 1^ and is due to the fact that stars enriched by a single 
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Fig. 3. a) A scatter plot of simulated stars in the [Mg/H]- 
[C/Mg] plane (Model I). The diagram shows 500 low-mass 
stars enriched by 1 — 20 SNe (Wi = 25). The yields are taken 
from WW95 and the masses of the SNe, ranging from 10 A4q 
to 100 Mq, are distributed according to a Salpeter IMF. b) A 
scatter plot of stars in the same plane as above, using yields 
from Netal97. It is seen that the existence of clear patterns in 
these A/H diagrams are dependent on the details of the pre- 
dicted SN yields 



SN have the lowest metallicities and the largest variations in the 
C/Mg ratio. When more and more SNe contribute to the metal 
content in the low-mass stars the metallicity increases and the 
variation in the ratio is averaged out. Note also that stars with 
a specific metallicity may have been enriched by quite a differ- 
ent number of SNe. For example, at [Mg/H]^ —3 the stars with 
lowest C/Mg ratio have been enriched by perhaps a couple of 
SNe while the ones with the highest ratio have been enriched 
by up to 17 SNe (see Fig. |a). 

Secondly, in Fig. ||a we see another effect. Instead of a 
decreasing scatter with metallicity the scatter is asymmetric, 
mimicking a trend. The C/Mg ratio seems to decrease with in- 
creasing Mg/H. This is not a normal evolutionary effect caused 
by time or metallicity (such as the decrease of [Mg/Fe] with 
[Fe/H] for [Fe/H]> —1 induced by the onset of thermonuclear 
supernovae (SNe type la) or the increase of [C/O] with [O/H] 
which could be explained by a metallicity dependent carbon 
yield as proposed by e.g. Gustafsson et al. 1999; Henry et al. 
2000) but rather a SN mass (i.e. number) effect. This is accom- 
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Fig. 4. a) Same as Fig. ^ except that the stars are coded. 
Circles (o) denote stars enriched by 1 — 4 SNe, squares (□) 
are enriched by 5 — 8 SNe, pluses (4-) by 9 — 13 SNe, dia- 
monds (O) by 14 — 17 and dots (•) are enriched by more than 
17 SNe. b) As above but with a different coding. Circles (o) de- 
note stars which have not been enriched by any SN with a mass 
above 27.5 Mq, squares (□) are enriched by one SN above 
27.5 Mq, pluses (+) have been enriched by two such SNe, di- 
amonds (O) by three and dots (•) by four or more SNe above 
27.5 Mq- The patterns displayed reflect the number of high 
mass SNe that have enriched the ISM before the sample star 
was formed 



plished by SNe producing a high C/Mg ratio at the same time 
produce a small amount of magnesium while SNe producing 
a low C/Mg ratio also produce much Mg. So, for an extremely 
metal-poor system an observed trend like this one does not nec- 
essarily imply chemical evolution in the normal sense (see also 
Tsujimoto & Shigeyama 1998). 

Thirdly, the stars in Fig. ^ja tend to group together in sub- 
structures. It is understood from Fig. ^ that these patterns are 
caused by a specific variation with progenitor mass in the car- 
bon yield. Roughly, one can say that these patterns are formed 
by a pronounced decrease in the yield around 27.5 Mq- We 
shall discuss these issues in more detail in Sect. ^ As we have 
mentioned, the substructures in Fig. |^a are sensitive to varia- 
tions in the mixing mass. As long as all star-forming regions 
have equal mass these substructures survive. This is unUkely, 
however. Nakasato & Shigeyama (2000) discuss metal enrich- 
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ment of the primordial ISM by individual SNe and find that 
the Mg/H ratio in various filaments may differ by ~ 1.0 dex, 
implying a large abundance scatter in the second generation of 
stars. Thus, a smearing effect most probably occurs in the hor- 
izontal direction in the diagrams of Figs. ^ and ^. The fact that 
we see patterns in the A/H diagrams is because we have not 
included this type of mixing in the models. 



3.4.2. A/ A diagrams 

Using Model I, it is possible to generate pure abundance ra- 
tios which can be displayed in A/ A diagrams such as those in 
Fig. ^. The patterns that are formed in these kinds of diagrams, 
opposed to the ones in the A/H diagrams, are insensitive to in- 
trinsic uncertainties such as mixing. They also show larger vari- 
ations in their shapes. For observational uncertainties of ~ 0.1 
dex in the Mg/H and C/Mg ratios the two scatter plots in Fig. ||a 
and b would appear quite similar apart from the different loci. 
However, for the same uncertainties in the C/Mg and Mg/Fe 
ratios the difference between the chemical patterns in Figs. ^ 
and |]d will survive. It is obvious that the stellar yields play a 
crucial role in the formation of these patterns. The other pairs in 
Fig. H show even larger differences due to larger disagreements 
in the two sets of yields (i.e. WW95 and Netal97). 

Model I and Model II generate A/A diagrams showing 
strong similarities. By comparing Fig. ^ and Fig. |^ we see 
that the stars are arranged in the same manner, apart from minor 
differences in the distribution within the formed pattern. This 
is a result of the discrete enrichment as the abundance ratios 
are not considered to be weighted with an IMF or some other 
mass distribution function. A set of yields generates a unique 
pattern while the mass frequency of SNe determines how the 
pattern will be populated by low-mass stars. In Sect. 4.2, we 
shall derive a mathematical expression for this statement. Note 
that the concentration of stars in the upper part of the distribu- 
tion in Fig. 1^ results from a higher density of individual SNe 
producing these abundance ratios (cf. the red curve in Fig. He). 



In Model II, the relaxation of the instantaneous recycling 
approximation within the star-forming regions introduces an 
evolutionary effect which has no counterpart in Model I. This 
effect should appear as a finite, non-vanishing dispersion in the 
different abundance ratios when the number of polluting SNe 
becomes large. The dispersion survives because stars that are 
formed early could only have been enriched by the most mas- 
sive stars while stars that are formed late (i.e. ^ 30 Myrs) could 
have been be enriched by SNe of any mass. However, the ef- 
fect can not easily be detected in the the patterns in Fig. |[ In 
fact, this is as expected since we only use a maximum number 
of 20 SNe in our simulations, which implies that the statistical 
fluctuations in the abundance ratios is still comparable to this 
dispersion. Eventually, the interaction between different clouds 
mixes the gas and the chemical abundance scatter in the subse- 
quent generation of stars is decreased. 



3.4.3. Possible sources of contamination 

In our models we have assumed that the only stars that en- 
rich a star-forming region with heavy elements have masses 
> 10 A^o, excluding the very massive stars (> 100 A4q). 
This is consistent with the life-time of the star-forming re- 
gions adopted, i.e. 30 Myrs. One could ask whether also 
intermediate-mass stars could be able to enrich a cloud. This 
would require that that either the life-time of the region is 
longer than assumed, or that the intermediate-mass stars were 
formed in epochs before the formation of the star-forming re- 
gion. In the latter case it is reasonable to assume that the region 
has also been polluted by more short-lived, high-mass SNe. 
That is, none of the Halo stars would be enriched by a single, 
intermediate-mass star The chemical patterns would be blurred 
but probably not very significantly. 

The patterns could also possibly be affected by a popula- 
tion of very massive stars, preceding the onset of normal core 
collapse SNe. Similarly to the intermediate-mass stars, the very 
massive stars would pollute the star-forming regions with un- 
known amounts of elements, not accounted for in the models. 

Another, related source to noise affecting the patterns is the 
possible existence of stars sampling the nucleosynthetic signa- 
ture of thermonuclear supernovae, like SNe type la. They are 
a different type of objects, not parametrized by the progeni- 
tor mass, and they will introduce a pattern in the A/ diagrams 
which is different from that of the core collapse supernovae. 
Progenitors of early thermonuclear SNe are thought to be close 
binary systems in which a white dwarf accrets matter from a 
subgiant star with a mass of 2-3.5 A4q (Branch 1998). The 
time delay between the onset of star formation and the forma- 
tion of these thermonuclear SNe are at least on the order of 0.1 
Gyrs. Thus, the core collapse SN patterns in the A/ diagrams 
may not be severely contaminated by these objects. 



As mentioned in Sect. 3.1 we do not consider metallicity 
dependent yields in our models. Such yields may also produce 
a smearing of the chemical patterns. The smearing is small for 
primary elements but could be as large as 0.5 dex for secondary 
elements, e.g., for the ratio i^N/^^'Mg (Umeda et al. 2000). 

In general, it would not be easy to recognize and remove 
stars enriched by these extra, hypothetical sources, except per- 
haps for stars with abundance patterns indicative of a pure ther- 
monuclear SN contribution. However, we have presently no 
strong reasons to believe that any of these sources have a sig- 
nificant effect on the A/ diagrams. 

3.4.4. Tine effects of statistical dependence 

Since the number of star-forming regions in reality is not in- 
finite there is a certain probability that two, randomly picked 
Halo stars may have formed in the same cloud, or equivalently, 
the observed sample of stars may not be completely statistically 
independent. If we randomly pick a number of stars which have 
been formed in a certain number of star-forming regions, how 
many of these regions are then represented by these stars? If 
the total mass of the Halo at the early epochs was ten times 
the mass of the stellar component today, which is ~ 10* Mq 
(Binney & Merrifield 1998), and all the gas was confined in 
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Fig. 5. Scatter plots of stars as they may appear in different A/ A diagrams, a-c) Simulated stars (Model I) enriched by 1 - 20 
SNe using yields from WW95. d-f) The corresponding A/ A diagrams produced with yields taken from Netal97. Note that for a 
particular abundance plane the difference (apart from the loci) between the distributions is large which arises from differences in 
the yields 



star-forming regions of 10^ Mq, there were approximately 
1000 different regions. Assume that every region formed an 
equal amount of stars. From available statistics (Christlieb & 
Beers 2000 and references therein) we estimate that there are 
approximately 50 000 extremely metal-poor Halo stars with a 
B magnitude < 17. Now, suppose that we randomly select a 
subsample of 100 stars. How many different star-forming re- 
gions (sampling different chemical series) are then represented 
by the sample? Given the much larger number of available Halo 
stars, we estimate that approximately 95% of the selected stars 
would sample different regions. On the other hand, if no more 
than 100 star-forming regions existed in the early Galaxy, about 
63% of the stars would still be statistically independent. 



We performed a small test to investigate the necessity of 
the assumption of statistical independence. We generated two 
samples containing 1000 stars each. The stars in the first sam- 
ple were selected from individual chemical series where no star 
was enriched by more than 25 SNe. These stars are statisti- 
cally independent. The second sample consisted of 40 chemi- 
cal tracks, i.e., all stars enriched by 1 — 25 SNe were selected 
from 40 different chemical series. Stars belonging to a chemi- 
cal track have a common chemical history and are statistically 
dependent. The difference between the two corresponding A/ 
diagrams was found to be relatively small and the characteris- 
tic pattern displayed by the first sample was well reproduced 
by the second one. Since only some tens independent chemical 
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Fig. 6. a) The distribution of stars in the [K/Sc]-[Sc/C] plane (cf. Fig. ge) as predicted by Model II. Due to the evolutionary effect 
the highly enriched stars do not converge to a point in the diagram even though they have been enriched by many SNe. This effect 
is not easily detected, however b) Same as Fig. |]: as predicted by Model II 



series contain enough information to form reliable patterns (re- 
member that 100 randomly selected Halo stars would sample 
maybe twice as many series), we conclude that the assumption 
of statistical independence is not vital for our conclusions. 

4. The origin of tlie chiemicai patterns 

In this section we shall discuss the origin of the chemical pat- 
terns found in the simulations above. We start with an approxi- 
mate approach to elucidate the source of the patterns, and also 
give some further insight into the reason for the variety of pat- 
terns, demonstrated in Fig. |] before we turn to the analytical 
theory. 

4. 1. Approximate approaches to pattern formation 

Let us, to be explicit, assume that the yields for three ele- 
ments, A, B and C, vary with mass such that, i.e, element A is 
predominantly produced by SNe in a certain progenitor mass 
range, mA> to niA", with yield pA', and for the rest of the mass 
interval the yield is assumed to be much lower and constant, 
PA" ■ Similarly, the element B is assumed to be produced by 
stars in the disjoint mass range to^' to m^", with yield ps', 
while the stars outside the interval produce the element B with 
the much lower yield pb"- The element C is assumed to be 
produced with a mass-independent yield pc- Denote the num- 
ber of SNe within the progenitor mass interval (mA','mA") by 
rii, and correspondingly n2 for the mass interval (mB',mB")- 
The total number of SNe is denoted na . 

With these assumptions it is easy to derive the main proper- 
ties of the distribution of stars in the [A/C]-[B/C] plane (Fig. 

We find the total yields of element A, B and C to be, re- 
spectively 
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Fig. 7. Generic A/ A diagram generated from artificial yields, 
two pronounced peaks in the A and B yields produce the char- 
acteristic "L"-shape with a cluster of stars round the centre of 
gravity 



{\og{pA" /pc)Aog{pB" /pc)}, and the end points of 
the "L" will be at {\og{pA' /pc)Aog{pB" /pc)} and 
{\og{pA"/pc)Aog{pB'/pc)}, respectively. These points 
correspond to cases with no SNe in the mass intervals 
(mA',rnA") and (mB'^mB"), (mB'^mB"), and (mA'^rnA")- 
Obviously, if these points can be observed, the yields can all 
be determined. 

If one further assumes that niPA' ^ (n-3 — and 
that n2PB' ^ ("-3 — "-2)pb" one finds the centre of gravity 
of the points in the A/ A diagram to be at {\og{pA' /pc) ~ 
log(n3/ni), log(pB//pc')-log(n3/n2)}. Thus, with the yields 
known observation of this point gives the ratios 
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For the standard deviation in [A/C] and [B/C] one finds 
with the assumptions above cr[A/c] ~ 0.4/^^777 and similarly 
<^[B/c] ~ OA/^. 

These estimates are, however, very dependent on the as- 
sumed yields - yields varying with stellar mass in more com- 
plex ways will contribute to the standard deviations, and thus 
make determination of ni and n2 in this way unrealistic. 

The two legs of the "L"-shaped distribution in Fig. ^result 
from the SNe with no strong contribution of elements A and B, 
respectively. The sparsely populated narrow sequences in Figs. 
H and ^ are, however, mainly the result of the pollution of the 
star-forming region by just one SN. The sequences are delin- 
eated by the range of SNe with different progenitor mass. This 
is illustrated in Fig. ^ where the corresponding distributions 
for just one SN (red dots), instead of 1 — 20 SNe as in Fig. ^ 
have been plotted. If such narrow sequences could be identi- 
fied observationally, which would require rich samples of very 
metal-poor stars and high observational accuracy, one might di- 
rectly read off the relative yields at different progenitor masses. 
The values of the latter will, however, remain unknown. 

In Fig. H we have also plotted the distributions of stars, re- 
sulting from cases with two polluting SNe (blue dots). It is clear 
from this figure, in comparison with Fig |[ that most of the 
structure of the A/ A diagrams is delineated already by models 
with two SNe. The density distribution in the diagrams of Fig 
||, is, however, determined by the IMF as transformed by the 
mass dependence of the yields. 



4.2. An analytical approach 

We shall now proceed to a more exact treatment of the forma- 
tion of the chemical patterns by deriving analytical expressions 
for the density distribution (i.e. frequency distribution) of stars 
in the A/ diagrams. This distribution is represented by a two- 
dimensional density function. It is constructed by a sum of den- 
sity functions, where each of these functions describes the dis- 
tribution of stars enriched by a certain number of SNe. In order 
to understand these functions and get a feeling for the param- 
eter dependence we shall begin by discussing one-dimensional 
density functions describing the distribution of a specific ele- 
ment or a ratio between two elements. The fundamental func- 
tions in this context are the ones that describe the distribution 
of stars enriched by individual SNe. In Appendix A we derive 
some general expressions for distributions of random variables 
and we shall frequently refer to those results in this section. 

The analytical theory is based on our simple model of 
chemical enrichment (i.e. Model I). Its only important parame- 
ters are the stellar yields and the IMF slope index. The chemi- 
cal patterns predicted by Model II are very much like the ones 
from Model I and the main result from Model I is not altered 
even if some changes in the density distribution of stars occur. 
The conclusion is that the stellar yields, or rather the variations 
in the yields with progenitor mass, play the crucial role for the 
shape of these density functions. 



4.2.1. The element distribution from individual 
supernovae 

We start with deriving the fundamental expression for the dis- 
tribution of stars enriched by individual SNe. In general terms, 
the relative frequency of stars in an A/ diagram is given by 
the relative frequency of heavy element producing phenomena 
as a function of the amount of heavy elements produced. In 
this study, we assume that these elements are produced in SN 
explosions and that the number density of SNe is basically de- 
termined by (j), the IMF. We have seen above that the actual 
number density is different from that of the IMF if, e.g., we 
consider continuous star-formation (Model II). However, this 
has not a big impact on the formation of the patterns. 

First, we need some definitions. The mass of a star formed 
in a star-forming region can be regarded as a random variable 
(r.v.) M, with a distribution function FM{m). The probability 
density function is then 



fM{m) = Flj{m) = (t>{m) (/)om 
normalized as 

■j){m)dm = 1. 



(10) 



(11) 



We say that M is an IMF(m,„in, mmax)-distributed random 
variable. Normally, the mass density, m(f){m), is normalized to 
one but since we are interested in number densities we choose 
to normalize 0(to) instead. Furthermore, let the yield be a con- 
tinuous function of stellar mass, x ~ p{m), defined on the 
interval [m„iin, m-max], i.e. [10, 100]. 

Now, as stellar masses are randomly distributed (according 
to this IMF) the amount of an element A, produced in a star 
of mass m, will also be randomly distributed {A stands for an 
arbitrary heavy element which is a product of stellar nucleosyn- 
thesis). However, the distribution of the r.v. ^(a) — P{A){M) 
is different from that of M. The probability of -'^(a) to be less 
or equal to x is gi ven b y the distribution function F(^a) (x), and 
according to Eq. ( A.l ) 



Fm{pI^^{x)), 



(12) 



assuming that p(A)(m) increases monotonically. For mono- 
tonically decreasing functions we will have that F(A)(a;) = 
1 — Fm{p'(^^j{x)). We shall write the element in the subscript 
within parentheses to distinguish these functions from the final 
expression which is written without parentheses. 

As our aim is to derive expressions for the density of stars in 
the A/ diagrams, we shall describe the distributions of random 
variables in terms of density functions. The density function, 
f(A){x), of X(A) is g iven by the derivative of F(^a){x) with 
respect to x. Eq. (A3) together with Eq. ( [lo| ) gives that 



X ^{pl2^{x)), (13) 



where p'^^-j = dpf^) / dm and p^^^) is the inverse function to 



(A) 



P(Ay Note that x e [X™", X^'^j'^], i.e. f(A)ix) as well as 
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Fig. 8. As Fig. |5]but for one and two polluting SNe only. The red dots denote stars enriched by one SN. These stars are located 
along one-dimensional curves in the abundance planes. The blue dots denote stars enriched by two SNe. Already two polluting 
SNe account for most of the structures present in the diagrams in Fig. |] 



P{A) {x) are defined on the interval between the minimum and 
the maximum produced amount of element A in stars with 
masses in the interval m e [m-min, m-max]- 

Eq. ( pj| ) holds for monotonic yields, P(a)- If the yield has 
local extre ma t here is no way of finding a single inverse to p(^) . 
As in Eq. ( |A.3[ ) it is necessary to split the interval of m such that 
the yield is monotonic on each subinterval. 

Let us consider a simple example. Assume that an element 
A is produced in massive stars in such a way that the stellar 
yield x = P(^)(to) oc to^/^ (measured in A4q). Furthermore, 
let the high-mass stars be distributed according to an IMF cx 
m~^. Thus, oc a;^ and oc m~^/^. Using 

Eq. ( [l3| ) we find that the number distribution of massive stars 
producing a specific amount of element A is proportional to 



(a;^) '^/{x^) ^1"^ ^ X ^. On a logarithmic scale the density 



function /(a) (2;) oc 10 according to Eq. (A. 16). 

The density function, ji^j^-^{x), describes the relative distri- 
bution of stars enriched by individual SNe, where each SN pro- 
duces a certain amount of each element. It consists of two fac- 
tors. The second factor, the IMF, accounts for the non-uniform 
distribution of stellar masses. It is a smooth, monotonically de- 
creasing function of mass, which, in practice, means that it will 
not produce any sudden changes in /(^j ix). On the other hand, 
the first factor, which depends on the stellar yield, may change 
drastically with mass. This is then reflected in the density func- 
tion. Fig. ^ shows two examples of density functions for the 
elements carbon and magnesium. The yields are given by the 
zero-metallicity models of WW95. 
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Fig. 9. a) The density function of carbon from individual SNe 
with the yield (Z = models) taken from WW95 (see Fig. |l]d). 
The Salpeter IMF is used. The histograms are binned data from 
numerical simulations of 20 000 clouds with a single SN explo- 
sion in each cloud (Model I) while the full line is the analytical 
functio n calc ulated from Eq. ( pj| ) and transformed according 
to Eq. ( A. 16 ). b) The corresponding density function for the 
element magnesium 



Note that we have not yet been discussing distributions of 
stars (i.e. low-mass stars) enriched in heavy elements. Thus, the 
variable x in /(y^) (x) is not the amount of an element in a low- 
mass star. So far, x represents the amount of the element pro- 
duced in the SN that has enriched the low-mass star In general, 
detailed knowledge about the mixing of the SN material with 
the ambient, possibly pre-enriched medium is needed in order 
to determine A/H ratios. In our simulations, this is accounted 
for by assuming a constant mass of the hydrogen clouds and 
no global mixing and we shall adopt the same mixing scenario 
here. This questionable assumption is of small significance for 
the pattern formation in the A/ A diagrams, i.e. for elements 
beyond hydrogen. 

4.2.2. The distribution from two or more supernovae 

If more than one SN explodes in each cloud the distribution 
of the freshly synthesized material in the low-mass stars is no 
longer described by the density function /(^)(a;). The sum of 
the contributions from every SN has to be considered, which 



Fig. 10. The two-dimensional yield ratio C/Mg as given by Eq. 
dl^. Also shown in the figure is the plane at P2(CMg) = 3.50. 
Every plane parallel to the mim2-plane which intersects the 
yield function at a certain point defines an integration region 
for which the integral in Eq. ( |T7[ ) can be calculated, giving the 
value of the distribution function in the point x = P2(CMg) 

alters the distribution. For example, two different sets of SNe 
may well produce the same total amount of an element. Thus, 
the correct density function describes a sum of the indepen- 
dent, equally distributed random variables, -'f (a) ■ This function 
is given by a convolution of /(y^) (x) with itself. 

The convolution formula for two random variables (i.e. two 
SNe) is derived in Appendix A. Here, we shall give the general 
expression for n random variables since we are interested in the 
shape of the density function when we sum up the contributions 
from n SNe. A generalization of Eq. (A. 7) is done by adding a 
third r.v. to the first two, then adding a fourth and so on. Thus, 



fn(A) (x) 



.f(A){xi)f(A){x2 - Xi) X ... 



X f{A)ixn-i - Xn-2)f{A)ix - a;„- 1 )da;i . . .da;^14) 

As for the expression in Eq. ( |A.7| ) the integrals are taken over 
the whole space. The function describes the distribution of the 
sum of n independent, equally distributed rv.s X^^j^y This is 
indicated by the notation n{A) in the subscript of fn(A){x)- 

It is worth mentioning a couple of properties of the con- 
volution described by Eq. (14). When n tends to infinity the 
corresponding density function of the arithmetical mean, i.e. 
the function derived from Eq. ( A.15| ), tends to a Dirac delta 
function centred at a: = / p{m)(j){m)dm. Moreover, a sum of 
random variables can be approximated by a normal distribution 
according to the Central Limit Theorem. The dispersion of the 
arithmetical mean is then proportional to n^^^^. 
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The number of SNe needed for this approximation to be 
valid depends ultimately on the shape of but n ~ 20 

is probably a good order-of-magnitude estimate. The averaging 
gradually erases the specific structures in the density functions 
and as soon as the shape of fn{A}{x) becomes well-behaved 
(i.e. Gaussian) it does not carry any information about the orig- 
inal yield and no new patterns are formed in the A/ diagrams 
(see, e.g., the decreasing scatter at the high-metallicity end in 
Fig. ^ or the crowding of stars round the point (—0.19, —0.25) 
in Fig. 

4.2.3. Abundance ratios 

In this section we shall derive the important expression for the 
distribution of an abundance ratio between two heavy elements, 
A and B. 

Let us first define a generalized stellar yield. The total 
amount of an element ejected by n SNe can be regarded as 
a yield in the n-dimensional m-space, i.e. 



p{mi 







(15) 



where p{m) is the normal, one-dimensional stellar yield for a 
star of mass m. In this way the yield ratio of two elements is 
simply defined as 



Pn{AB){'mi, ...,TO„) 



E P(B){ini) 



(16) 



As before, the subscript n stands for the number of polluting 
SNe which in this case is equal to the number of dimensions as 
the original stellar yield is one-dimensional. 

No w, si milarly to the multi-dimensional random variable 
in Eq. (^JJ), we form the r.v. Xr,(AB) Pn{AB) {Mi, il/„). 
This ratio has properties similar to that of the mean of a single 
element A. It is possible to show that the distribution function 
Fn{AB) {^) can be written as 



FniAB){x) 



(mi) X ... X 0(m„)dmi...dm„ (17) 



according to Eq. ( ). The integration region is defined by 
the equation 



Pn{AB){mi, ..^rrin) < X, 



(18) 



i.e. the part of the mi...TO„ -hyper plane where the yield is less 
than or equal to x (see Fig. HQ). Note that because all stellar 



mass r.v.s are independent the integrand in Eq. ( A. 9 ) reduces to 
a simple multiplication of the IMFs. 

The integrand in Eq. (|l7|) is almost trivial to handle while 
the integration region is highly non-trivial. From Eq. (^s]) we 
see that depending on how the yield varies, may not even 
be simply connected. This is illustrated in Fig. The inte- 
gral should be calculated for every x which then gives the dis- 
tribution function. The density function is obtained taking the 
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Fig. 11. a) The density function of low mass-stars in [C/Mg]. 
Each star has been enriched by 3 SNe. Their masses were 
distributed according to the Salpeter IMF. As before, the his- 
togram shows binned data from a numerical simulation while 
the full line is the corresponding density function to the formal 
solution of the distribution function in Eq. (p^). b) The density 
function for [Mg/H]. This is in principle a convolution of the 
random variable ^(Mg) (see Fig. ^) using Eq. (jl^) for n = 3. 
The magnesium abundance is measured relative to 5 x 10^ -Mq 
of hydrogen 



derivative with respect to x as in Eq. ([l3|). Fig. ITT] shows two 
examples of density functions of such compound random vari- 
ables. Fig. [ri| a shows the distribution of stars over [C/Mg] for 
three polluting SNe using the same yields as above. The full 
line is a solution to Eq . (p^ . The density function in Fig. pl] b 
is calculated from Eq. (|14j). Both functions are transformed to a 
logarithmic scale (and translated relative to solar values) using 
Eq. (O^. 



4.2.4. The construction of A/ diagrams 

We are ultimately interested in two-dimensional distributions 
describing the density of extremely metal-poor low-mass stars 
in the A/ diagrams. If two random variables are independent 
their joint density equals the product of the individual densities 
(see the derivation of the convolution formula in Appendix A 
where independence is assumed). However, the random vari- 
ables we discuss here are dependent. If we simultaneously ob- 
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Fig. 12. Intersecting integration regions. ADy (dark-shaded 
areas) is defined by P2(CMg) — V ^ 3.50 (cf. Fig. [l^ 
while ADx and ADx' (light-shaded areas) are defined by 
P2(Mg) = X ^ 0.10 and P2(Mg) — x' ^ 0.20 respec- 
tively. Hence, AD^y and AD^'y are given by the intersections 
AZJj; Pi ADy and AD^' H ADy. A calculation of the integral 
in Eq. (|9|) over AD^y for every a; = P2(Mg)('™ii '^12) and 
y = ?'2(CMg) ("^i7 "^2) gives the two-dimensional density func- 
tion /2(Mg,CMg) 



serve the random variables in Fig. |1 l| a and b the corresponding 
density function is not a multiplication of the two individual 
density functions. It is not possible to obtain every value of the 
ratio [C/Mg] for a given value of [Mg/H] by the combination of 
three SNe as the two variables are entangled via the progenitor 
masses of these SNe. 

Suppose that we have four elements A, B, C, and D, where 
not all elements are necessarily different. Now, the joint distri- 
bution of X„(^B) and X^i^cd) (or, e.g., and X„(bc')) is 
described by the function fn(AB,CD) (oJ'fniA,BC))- As usual, 
the subscript n stands for the number of contributing SNe. The 
derivation is similar to that of the one-dimensional distribution 
function in Eq. ( [l7|). He re, the density functions are calculated 
directly, as in Eq. (|A.14). Thus, 



In{AB,CD){x,y) = / ■■• / (/)(mi) X ... 

X 0(m„)dmi...dm„, (19) 
where the integration region, AD^y is defined via the equations 



AD^: X <Pn(AB){mi,...,mn) <x + Ax 
ADy-. y <Pn^cD){mi,...,mn) <y + dy. 



(20) 



Thus, for every point [x, y) the integration is performed over 
the small common region in m-space defined by the intersec- 
tion, ADx n ADy. The density of stars (enriched by n SNe) 
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Fig. 13. Contour plot of /Mg,CMg (cf. Fig. ^) computed for 1- 
20 SNe with w, = 1/20 = 0.05 using the yields from WW95. 
The IMF is of Salpeter form. The relative number density of 
stars is indicated by the colour bar 
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Fig. 14. Contour plot of /ksc,psc- The model parameters are 
the same as in Fig. |l3[ The stellar yields are taken from 
Netal97. The colour bar indicates the relative number density 
of stars 

in the [A/H]-[i?/C] plane, described by fn(A,BC) {^i y)^ is cal- 
culated similarly. However, the integration region is changed to 
(see Fig. |12|) 



AD.^ : X < Pn{A){mi, ... 
ADy : y < Pn{BC){fni, ■ 



m„) < X + da; 
■,mn) <y + dy. 



(21) 



The integrands in Eq. ( ]17[ ) and Eq. (g_9|) are the same. The 
difference lies again in the integration regions and we note that 
the (generalized) stellar yields are important components in the 
expression for the abundance ratio as the functions defining 
these regions. Moreover, the IMF only partly determines the 
density in each point, together with the gradient of the yield 
function. This is clearly seen in Eq. ( [l3| ) but can also be traced 
in Fig. |2| where the differential areas have different sizes de- 
pending on the shape of the yields. 

The final step to an analytical expression for the distribu- 
tion of stars in an ^/ diagram is a summation of the density 
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functions over the number of SNe. The expression for the den- 
sity in an Aj A diagram, i.e. the density in the \AIE\-\CID\ 
plane, is 



fAB,CD{x,y) X fi{AB,CD){x,y). 



4=1 



The weight Wi is the fraction of clouds forming i SNe. If 
we substitute fi(A,BC) v) ™ the sum, we instead obtain a 
function formally describing the density of stars in the [A/H]- 
[B/C] plane. Note that the first term in the sum is a one- 
dimensional function since any star enriched by a single SN 
has a unique chemical composition determined by the yields of 
that SN. Thus, the density function maps the one-dimensional 
yield ratio on the A/ A plane. The function is seen as a curve in 
Figs. pj|and |l^(see also Fig. The form of this curve actu- 
ally determines the form of the total density function to a large 
extent. 

We end this section by a small note on scale transforma- 
tions. The appropriate equations are derived in Appendix A. In 
our derivations all functions are on a linear scale. Often astro- 
physical quantities may span several orders of magnitude. It is 
then more advantageous to display these quantities on a log- 
arithmic scale. The transformation is made using Eq. (A^.16) 
and Eq. (A. 17) for one- and two-dimensional functions respec- 
tively. However, it is also possible to directly take the logarithm 
of the yield or the generalized yield and use these functions in 
the expressions for the densities instead. 

4.3. The robustness of the A/ A diagrams 

It is not probable that patterns like those in Fig. |l3| can be 
observed, since they are based on oversimplified models. The 
star-forming regions have different masses. Mixing within the 
regions is not complete and global mixing of remnant, enriched 
gas, or infalling gas, and the formation of a second genera- 
tion of star-forming regions probably occurred even for the ex- 
tremely metal-poor Halo. By neglecting all these effects the 
treatment of the mixing in our models is oversimplified (see 
Nakasato & Shigeyama 2000). Thus, a horizontal smearing ef- 
fect occurs in A/H diagrams such as [Mg/H]-[C/Mg], and the 
patterns are most likely lost in these diagrams. 

The patterns in the A/ A diagrams are, however, not sensi- 
tive to the masses of the star-forming regions, i.e. the amount 
of ambient gas that is mixed with the SN ejecta. If the SN rem- 
nant material is mixed with a pre-enriched ISM one might yet 
think that the information on the original production sites of 
the elements can not be traced in the observed abundance ra- 
tios, as different fractions of the enriched medium may con- 
tribute differently to different stars. This occurs when a second 
generation of star-forming regions is formed. Suppose for ex- 
ample that a new region is formed out of the dispersed gas of 
two former such regions with the proportions 1 : 3, i.e. 25% 
of the total mass originates from the first region and 75% of 
the mass from the second region. Assume further that these 
former regions were enriched by one SN each and the new re- 
gion produces one SN. An abundance ratio in the gas and in 
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Fig. 15. Contour plots of the density of stars in the [Na/Mn]- 
[Mn/Co] plane with yields from Netal97. a) The density func- 
tion with unweighted yield terms (i.e. ~ 1 for all i). b) 
The corresponding function with weighted yield terms accord- 
ing to Eq. ( p3[ ) representing the effect of global mixing. The 
weights have a uniform distribution, ranging from to 1. Due 
to the weights the microstructure in the second A/ A diagram 
has been wiped out. Otherwise, the difference between the two 
diagrams is very small 



subsequently formed stars is then found by summing up the in- 
dividual contributions from each SN, as in Eq. (|l6|). However, 
each term has a a weight associated with it, depending on how 
great the contribution is from each region. In our example the 
weights are 0.25, 0.75 and 1.00. Generally, Eq. ([l|) has to be 
modified such that 



E i>.P(B){mi) 

i=l 



(23) 



The only difference between this expression and the one in Eq. 
^ is the weight on each yield. The weight is supposed to 
mimic the effect of large-scale mixing and turbulent motions in 
the ISM and can be regarded as a random variable similar to 



the mass such that the rv. X 



n(AB) 



IS 



X 



n(AB) 



E 2iJ3(^)(M,) 
i=l 

n 

E ^tP{B){M^) 



-p«(^^)(Si,Mi,...,S„,M„)(24) 



where the superscript ^ denotes the fractional contribu- 
tion defined by the weighted generalized yields in Eq. 
(p3|). By introducing these weights, the generalized yield 
Pri(A_B)(™i- ■■■,'mn) looscs its Symmetry (e.g. in Fig. the 
two-dimensional yield ratio would no longer be symmetric 
around mi = 7712). Thus, depending on the weights, a specific 
abundance ratio may be produced by several different sets of 
SNe. This means that if we have no information on the weights 
we are no longer able to identify the set of SNe that gives the 
observed ratio. 
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Even so, simulations with randomly distributed weights 
show that the distribution of stars in A/ A diagrams is practi- 
cally unaltered (see Fig. |l^ with respect to the original distri- 
bution for which all = 1. This is due to the fact that the in- 
tersection between the two regions AD| and AZ?^ is non-zero 
(thus, the integral in Eq. ( p^ ) is non-zero) in almost the same re- 
gion in yield space (i.e. abundance ratio space) as for the unper- 
turbed (unweighted) integration region ADxy Note however, 
that the integration region over stellar masses may be different 
and far from symmetric which changes the value of the density 
function in each point. This is also observed in Fig. |l5|. The 
distributions have the same shape but the density of stars in the 
structures is somewhat different. Furthermore, in Fig. |ljb the 
microstructure is wiped out. It is clear that some information 
must be lost by introducing weights on the yield terms. Thus, 
we have no longer knowledge of the absolute amounts of the 
elements produced in the SNe. 

Even though the two chemical patterns in Fig. [ij look qual- 
itatively similar a relevant question is whether the microstruc- 
ture in Fig. [l^ a is crucial for a reconstruction of SN yields. This 
structure originates from the one-dimensional curve sampling 
individual SNe (cf. the red curves in the diagrams of Fig. ^) or 
sooner it corresponds to this curve for the case n ~ 2,3, etc. In 
principle, this structure does not contain any unique informa- 
tion and the loss of it should not be vital for a yield reconstruc- 
tion procedure. In this particular example the pattern in Fig. 
p3| a also shows two distinct peaks of which the right one is the 
convergence point or the centre of gravity. This dichotomy is 
not as prominent in Fig. Would the lack of a second peak 
lead to a significantly different result when reconstructing the 
yields? We believe that this is not the case. A comparison with 
the much greater diversity of patterns shown in Fig. ^ suggests 
that the major part of the information content on the yields in 
Fig. [l5| a still remains in Fig. [l5| b. A number of similar simu- 
lations for other elements confirm this result. However, more 
tests should be carried out in order to quantitatively determine 
the sensitivity of the yields on the details of the chemical pat- 
terns. 

5. Discussion - Can patterns be observed? 

Above, we have discussed the formation of patterns in A/ dia- 
grams and the possibilities of using abundance ratios as a diag- 
nostic of early stellar nucleosynthesis. However, the question 
remains: Are chemical abundance patterns detectable in prac- 
tice? 

It is difficult to estimate the accuracy needed in the abun- 
dance determinations and the minimum number or sample stars 
that is required for detecting the chemical patterns. Since the 
number of Halo stars that can be observed accurately is lim- 
ited, these parameters are dependent. Fig. 16 shows a compar- 
ison between the original distribution of about 100 stars in the 
[K/Sc]-[P/Sc] plane (cf. Fig. |l^ and the corresponding distri- 
bution for which an uncertainty of 0.08 dex is added in the 
observed abundance ratios. The original "V"-shape is barely 
detected in Fig. |l6| b. However, the pattern would more easily 
appear if we could double the number of stars. Larger uncer- 
tainties in the abundance ratios would completely wipe out the 




2.0 



1.6 



1.4^ 
0.2 



(b) 



0.4 



0.6 
[K/Sc] 



0.8 



1.0 



Fig. 16. a) A scatter plot of around 100 stars (enriched by 1-15 
SNe) in the [K/Sc]-[P/SC] plane. The nucleosynthesis calcula- 
tions by Netal97 have been used. The "V"-shape is easily de- 
tected (cf. Fig. |l^. b) The corresponding distribution with an 
additional uncertainty of 0.08 dex in the abundance ratios. The 
"V"-shape is less pronounced but still detectable 



pattern and this can not be compensated for by adding more 
stars. On the other hand, the "U"-shape in Fig. ||c is detectable 
even for uncertainties above 0.1 dex in the abundance ratios. 
However, for making that pattern visible at all several hun- 
dred stars must be observed, many more than are needed for 
the "V"-shaped pattern. Thus, depending on the shape, some 
patterns are more sensitive to observational uncertainties than 
to the number of sample stars while the opposite is true for 
other types of patterns. 

This discussion also holds for A/H diagrams (e.g. Fig. |]). 
Assuming that there are no intrinsic uncertainties which wipe 
out the patterns, we estimate, by convolving the density func- 
tion in Fig. [l^ with a Gaussian profile, that the structures be- 
come undetectable if the uncertainties in the derived abun- 
dances are larger than ^ 0.05 dex. This is an upper limit as 
we in reality only deal with a limited number of sample stars, 
i.e. poorer statistics. However, e.g., the bimodal structure in the 
[N/0]-[0/H] plane (see Fig. 3 in Karisson & Gustafsson 2000) 
caused by a pronounced peak in the N-yield (WW95) is not 
erased even for uncertainties in the abundances as large as 0.15 
dex. This also suggests that large-scale structures in A/H dia- 
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grams can possibly survive the effect of mixing in the interstel- 
lar medium, if the mass range of typical star-formation regions 
is Umited to within a factor of 2 or so. One should also note that 
similar patterns could also be expected to be visible in smaller 
star systems enriched by a finite number of SNe with full mix- 
ing of the system occurring between each SN. 

Observing chemical patterns will be a challenging task. 
With the tools of today and current methods for abundance 
analysis we are able to decrease the absolute observational un- 
certainties in the abundance ratios to ~ 0.1 dex. A strictly dif- 
ferential study may reach the 0.05 dex level of uncertainty. 
This may be slightly too large to allow detection of the fine- 
structures in many patterns which seem to begin appearing first 
at a level of about 0.05 dex. However, larger uncertainties can, 
as we have seen, partly be compensated for by a larger stellar 
sample. Also, in a number of cases patterns should be visible 
already with errors in the relative abundances of ^ 0.1 dex. 

5. 1. Existing evidence for ctiemical patterns 

Our aim is here not to make a literature survey of abundance 
data for metal-poor stars or to re-analyse existing data but 
merely to point out some studies and phenomena that possible 
can be related to variations in stellar yields and the formation 
of chemical patterns. 

There may be undetected patterns in the sample of some 
thirty metal-poor giants ([Fe/H] < — 2) observed by McWilliam 
et al. (1995). The sample is fairly homogeneous although the 
stars are not dwarfs and internal nucleosynthesis may have al- 
tered some of their surface abundances. Nevertheless, let us 
look at the distribution of stars in the [Mg/Fe]-[C/Mg] plane 
which is shown in Fig. |l^ The symbols are shaded according 
to the stellar metallicity as measured by [Fe/H]. The dispersion 
is large in both directions, spanning 1.5 dex in [Mg/Fe] and 
nearly 3 dex in [C/Mg]. The distribution of stars is asymmet- 
ric and there is no strong correlation with metallicity. However, 
there are too few stars in the sample to really allow the detec- 
tion of any patterns. 

If we instead plot the stars in the [C/Ca]-[C/Mg] plane, we 
find something completely different. The stars form a quite 
beautiful relation (see Fig. 18). The slope is close to unity 
and there seems to be no direct dependence on the iron abun- 
dance as indicated by the shading of the star symbols. In accor- 
dance with our discussion on early chemical enrichment we 
shall assume that these stars have been enriched by a small 
number of core collapse supernovae at the epoch of forma- 
tion of the Galaxy. If so, this leads us to beheve that, what- 
ever the variation of the yields with SN mass is, the ratio of the 
Ca-yield to the Mg-yield is rather independent of progenitor 
mass. The scaling factor (by number relative to solar values) 
can directly be estimated from the offset in Fig. [T^ and gives 
Nca — l.lSiVMg although it is consistent with unity. This re- 
lation can also be observed in the [Fe/H]-[Ca/Mg] plane where 
the stars are scattered around the constant value of 0.06. The 
scaling relation between the yields (by mass) is then estimated 
topca(w) ~ O.lpMg (to)- Only weak observational constraints 
on the carbon yield can be deduced from the diagram in Fig. 
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Fig. 17. Observed abundances of stars from McWilliam et al. 
(1995) displayed in the [Mg/Fe]-[C/Mg] plane (cf. Figs. ^ and 
||d). The shading according to the bar on the right indicates the 
metallicity measured by [Fe/H]. Stars with B — V > 0.80 are 
removed from the sample in order to minimize the effect of 
carbon depletion in the red giants. There are too few stars to 
allow detection of any patterns although there seems to be a 
large, asymmetric dispersion 



|18[ However, the observed abundance ratios span almost three 
orders of magnitude which is considerable. 

Recently, Jehin et al. (1999) observed a sample of mildly 
metal-poor stars and found interesting correlations between the 
abundances of the a-elements and the r- and s-process ele- 
ments. The a-elements were correlated with the r-process el- 
ements in a one-to-one relation while a subpopulation of the 
stars with high [a/Fe] ratio seemed to be enriched in the s- 
process elements (see their Fig. 7). This excess in s-process 
elements was interpreted as an accretion phenomenon in dense 
environments, presumably globular clusters. We note that their 
correlation diagrams show strong similarities with some of the 
A/ A diagrams presented here, see Fig. ^ and Fig. ^ There is 
little doubt that the reversed "L"-shape and their two-branches- 
pattern are both caused by strong variations in the stellar yields. 
However, it is not clear whether the two-branches-pattern could 
be explained by yield variations alone or if an evolutionary ef- 
fect is necessary as proposed by these authors. 

Another possibly related phenomenon is the existence of 
CN-strong and CN-weak line stars in some globular clusters 
(see e.g. Cannon et al. 1998). There seems to be a distinct bi- 
modal abundance pattern in this population of stars, not much 
different from patterns formed by a yield with a pronounced 
maximum. The bimodality in the CN line strengths is not yet 
understood but, again, strong variations in the yields would pro- 
duce similar patterns. We should emphasize that our discussion 
on the formation of patterns only holds in a strict sense for ex- 
tremely metal-poor environments while the latter two examples 
concern relatively metal-rich systems. 

6. Conclusions 

Observations of the most metal-poor Galactic halo stars show 
convincing evidence for a large star-to-star scatter in abun- 
dances relative to hydrogen as well as abundance ratios for a 
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Fig. 18. Abundances of stars from the sample of McWilliam 
et al. (1995) plotted in the [C/Ca]-[C/Mg] plane. As for Fig. 
[l7[ , the symbols are shaded according to the stellar metallicity. 
Stars with B — V > 0.80 have been removed from the sample. 
Note that the observed stars form a tight relation with a slope 
close to one 



variety of elements, not the least the neutron-capture elements 
(McWilHam et al. 1995; Ryan et al. 1996; McWilliam 1998; 
Burris et al. 2000). This phenomenon can most easily be ex- 
plained in terms of local enrichment of the primordial ISM by 
a small number of exploding massive stars (Audouze & Silk 
1995). Since the amounts of newly synthesized elements de- 
pend strongly on the mass of the exploding star (some elements 
are also affected by parameters such as rotation and metallic- 
ity), the abundances might have varied extensively through- 
out the Halo ISM. Before turbulent motions in the ISM had 
time to wipe out the chemical inhomogeneities, formation of 
low-mass stars occurred and the inhomogeneities could be pre- 
served. Hence, studying these stars by statistical means, espe- 
cially by displaying them in A/ A diagrams relating different 
abundance ratios, reveals important information on the produc- 
tion sites of the elements. 

We have demonstrated that a sample of extremely metal- 
poor stars displayed in A/ diagrams forms patterns which orig- 
inate from specific variations in the stellar yields. Thus, the 
formation of patterns is a natural consequence of the varia- 
tions in the amount of SN-produced material. This is clearly 
seen in the analytical theory. The form of the patterns depend 
critically on the shape of the integration regions which are de- 
fined by the SN yields. Furthermore, assuming that the ejected 
matter from SN explosions can be considered chemically ho- 
mogeneous (see e.g. Kifonidis et al. 2000; see also the com- 
ment by Arnett 1999) we claim, based on the result given in 
Sect. 4.3, that chemical abundance patterns in A/ A diagrams 
(such as [C/Mg] vs. [Mg/Fe]) survive the effects of large-scale 
mixing in the ISM (see Fig. |5]). Moreover, a comparison be- 
tween the simulations of Model I and Model II (cf. Fig. |^ and 
Fig. ^ indicates that the patterns are not very sensitive to the 
mass distribution function of exploding SNe (see also Fig. 3 
in Karlsson & Gustafsson 2000). In the analytical expression 
the mass distribution function appears in the integrand. Hence, 
it merely specifies the relative density of stars within the pat- 
terns, not the form of the patterns. For a finite stellar sample. 



however, the apparent form of the patterns will depend slightly 
on the mass frequency of SNe since a finite number of stars do 
not cover all possible values defined by the theoretical density 
function. 

Continuous star formation within the clouds (Model II) 
could modify the chemical patterns. However, such an evolu- 
tionary effect would result in a displacement of the centre of 
gravity rather than the formation of new patterns. The number 
of contributing SNe in individual star-forming regions must be 
large (> 20) in order to detect the displacement (cf. Figs. 
^ and Fig. ^ and it vanishes for stars formed in a second gen- 
eration of star-forming regions. The general characteristics of 
the patterns are, however, not affected (see Sect. 4.3). 

In comparison to studies of individual Halo stars with dra- 
matic abundance signatures hke CS 22892-052, the inclusion 
of less extreme stars makes it possible to quantitatively discrim- 
inate between chemical patterns formed by different SN yields 
(see Fig. Future prospects include studies of large, homoge- 
neous samples of extremely metal-poor Halo stars with accu- 
rately determined abundance ratios. In fact, several surveys of 
this kind are already under way. The idea would be to recon- 
struct the yields by solving the inverse problem. An observed 
chemical pattern in some A/ A plane is compared with simula- 
tions whereupon the yields are changed until the two distribu- 
tions are, in some sense, equal. Some potential complications 
are noted. The possible contamination of thermonuclear SNe 
(SNe type la), intermediate-mass stars, and/or very massive 
stars might affect the chemical patterns and blur the nucleosyn- 
thetic signature of the core collapse SNe. Metallicity dependent 
yields affect the patterns of secondary elements. However, we 
estimate that for the first generations of extreme Pop. II stars, 
and for several elements, the contamination should not be very 
severe. The development of more realistic, time-dependent 
models should elucidate these problems, as well as those of 
mixing on various spatial and temporal scales. In a forthcom- 
ing paper we shall present a general method for SN-yield re- 
construction. 

A quantitative analysis is possible only if the SN yields can 
be strictly parametrized. In particular, the explosive nucleosyn- 
thesis is a product of a chaotic behaviour of the stellar material, 
intimately connected to hydrodynamical instabilities. Shocks, 
convection and turbulent motions in the SN gas could introduce 
intrinsic uncertainties in the produced amount of elements. The 
stochastic nature of such an effect could blur the chemical pat- 
terns in the A/ diagrams. However, if the effect is sufficiently 
small it should be possible to retrieve mean yields as functions 
of mass, rotation etc. We conclude that the chemical patterns 
are useful diagnostics of yields of core collapse SNe, and if 
patterns would be detected, we would have the opportunity to 
probe the earliest phases of stellar nucleosynthesis. 
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Appendix A: 



If Y and Z are independent Eq. ( 4.4 ) reduces to 



The analytical treatment of the stellar distributions in the A/ 
diagrams (Sect. ^ is developed in order to form a deeper un- 
derstanding of the parameter dependence of the observed pat- 
terns. Here, we shall derive some general expressions for den- 
sity functions as well as distribution functions corresponding 
to an arbitrary random variable or pairs of random variables. 
Note, that in Sect. 4.2 we wrote the random variables in the 
subscripts within parentheses. This notation was adopted for 
purposes explained in the text. However, this in not the con- 
ventional notation, and in this section the subscripts are written 
without parentheses. 

A.1. Functions of one random variable 

Suppose we have a random variable (r.v.) X and a function 
y — g{x). Assume this function to be monotonically increasing 
with increasing x. Now, define the rv. Y ?&Y — g{X). The 
probability that Y < y'K then given by the distribution function 
forF, 



Fx{x) 



fY{v)fz{.z)dydz. 



(A.5) 



y+z<x 



This expression can be rewritten as 



+ 00 



Fx{x) 



fviy) 



fz{z)dz 



dy 



+ 00 



fY{y)Fz{x - y)dy. 



(A.6) 



The derivative of Fx{x) with respect to x gives the density 
function. 



fx{x) 



fY{y)fz{x - y)dy 



(A.7) 



FY{y) = P{Y <y) = P{g{X) < y) = 

= P{X<g-\y)) = Fx{9-\y)), 



(A.l) 



where Fx is the distribution function of the r.v. X. The inverse 
function of g is denoted by g^^ . For a monotonically decreas- 
ing function g we will have that Fyiy) = 1 — Fx{g^^{y)) 
instead. 

In the continuous case the corresponding density function 
of Y is given by the derivative of Fyiy) with respect to y such 
that 



fYiy)^F^ig-\y))^fxig-Hy)) 



dy 



(A.2) 



We can allow the function g to be non-monotonic by defin- 
ing it as a sum, g — gi + ...+gn where each function gi is mono- 
tonic and equivalent to g on the open subinterval ]xi-i,Xi[ and 
zero elsewhere. The x^, z = 1, n — 1 are real roots of g 'jx) 
and xq, Xn are the end-points. Thus, similarly to Eq. (A.2) the 
expression for fviy) reads 



(assuming that derivation under the integral sign is allowed). 
Eq. (A^) is the convolution formula for two continuous and 
independent random variables. 

A.3. One- and two-dimensional distributions ofn 
random variables 

In this section, we shall form expressions for one-dimensional 
density functions of n continuous random variables as well as 
two-dimensional density functions of n random variables. 

Given n random variables Xi,...,Xn and a real-valued 
function g — g{xi, Xn) we may form the one-dimensional 
r.v. 



X = 9(^1: ■■■,^n)- 



(A.8) 



For every given number x, we denote by the region in the 
xi...Xn hyper-plane such that g{xi, ...,a;„) < x. Hence, the 
distribution function of X is given by the integral 



fYiy) = :r-FY{y)-Y: 



{97\y)) 



fx{9-\y)) 



dy ^ dy 

hHis-\y))y ^ ^ ^ 

where g' = dg/dx. 



A.2. The convolution formula 

Given two continuous random variables, Y and Z, with distri- 
bution functions Fy and Fz respectively, the sum, X — Y + Z, 
is described by the distribution function Fx given by the inte- 
gral, 



Fxix) 



'y+z<x 



lY,z{y,z)dydz. 



(A.4) 



Fx{x) = P{X < x) = P{{X^, ...,X„) e D,) = 

= ■■■ fxu.-,x^{xi, ■■■,Xn)dxi...dXn, (A.9) 

J Jd^ 

where fx^ x„ {xi, ■■■,Xn) is the density function of the joint 

distribution of Xi, X„. The corresponding density function 
of the r.v. X is simply give n by the derivative of the distribution 
function Fx{x), a s in E q. (A.2). 

Clearly, Eq. (A.9) is in particular valid in the one- 



dimensional case. Namely, if t hat y = g{x) is a monotonically 
increasing function, then Eq. ( A^) reduces to Eq. ( [4.l| ). 

Again, suppose that we have n random variables. It is then 
possible to define two new r.v.s X and Y via the functions 

g{xi, Xn) and h{xi, Xn) such that 



X 



Y = h{Xi,...,Xn). 



(A. 10) 
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The expression for the joint distribution of X and Y is simi- 
lar to that of Eq. ( |A.9| ). However, we shall directly express the 
density function fx.vix, y) by defining the integration region 
as an intersection between the two differential integration re- 
gions ADx and AZJj,, 



{x<X<x + dx}^ {(Xi, X„) e Ai?J 
and 

{y<Y<y + Ay} = {(Xi, X„) e ADy}. 



(A.ll) 



(A. 12) 



Note, that the density function of X, defined by Eq. ( |A.8[ ), can 
be calculate d dir ectly by substituting AZ)^; for in the inte- 
gral in Eq. ( |A.9| ). Now, for calculating the joint statistics of X 
and Y we introduce AD^y — AD^ f] Ai'y, or 



(A.13) 



{x < X < X + dx, y < Y < y + dy} = 
^{iXi,...,Xn)eAD,y} 

Finally, we then have 

fx.Y{x,y) = 

= P{x < X <x + dx,y <Y <y^ dy) = 

- P((Xi,...,X„) e AD„j) = 

= / •■■ / fxi,...,x^{xi, ■■■,Xn)dxi...doiAM) 

A general treatment of multi-variate distributions is dis- 
cussed in e.g. Papoulis (1991). Further reading on the theory 
of probability can also be found in Gut (1995). 

A.4. Some scaling laws 

In connection to Sect, i.2.2 and Eq. (|l^ we would also like to 



consider a density function describing the mean of rv.s instead 
of the sum. With a simple scaling /„x (x) is transformed into 



f{nX)ix) = nfnxinx), 



(A.15) 



where the subscript (nX) on / denotes the average of n ran- 
dom variables. 



All functions in Sects. 4.2 and 4.3 are displayed on a log- 
arithmic scale. Using Eq. (A.l) and Eq. (A.. 3) with g{X) — 
\og{X) we see that an arbitrary function of one variable trans- 
forms as 



/iogx(x)=ln(10)xlOVx(10^). 



(A. 16) 



The factor In(lO) arises from the fact that we use the 
10-logarithm, not the natural logarithm. Similarly, a two- 
dimensional function transforms as, 

/iogxaogy(a;,y) -ln2(10) x 10^+Vx,y(10", 10^). (A.17) 

However, instead of calculating the densities on a linear scale 
and then make the transformations, it is easier to directly take 
the logarithm of the generalized yields and use these functions 
in the calculations. 
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